Three dimensional topological invariants for time reversal invariant Hamiltonians and 

the three dimensional quantum spin Hall effect 
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The Z2 invariant for filled bands in the ground states of systems with time reversal invariance 
characterizes the number of stable pairs of edge states. Here we study the Z2 invariant using band 
touching methods discussed in a recent previous work Q and extend the study to three dimensions. 
Band collisions preserve the Z2 invariant both in two and three dimensions, but there are crucial 
differences in the two cases. In three dimensions, we find a novel fourth Z2 invariant which is 
characterized by a "trapped monopole" in momentum space. If the monopole charge in half the 
Brillouin zone is odd, then atleast one of the monopoles cannot recombine with another monopole 
and vanish unlike the case when the monopole charge is even. We also point out the possibility 
of a three dimensional quantum spin Hall effect and discuss the connection of various topological 
invariants to such an effect. 

PACS numbers: 



INTRODUCTION 

A new invariant that characterizes ground states of two 
dimensional systems in a periodic potential with time re- 
versal invariance when the Fermi Energy lies in a gap was 
introduced recently by Kane and Melc 2J. They discov- 
ered this invariant while studying a model for graphene 
which exhibits the quantum spin Hall effect |3j. The 
model has a single pair of edge states which are robust to 
small perturbations. The robust edge states are also seen 
in an earlier model for the quantum spin Hall effect based 
on Landau levels in the presence of a strain gradient |J| 
and a model for the square lattice [5J- The robustness 
of an odd number of time reversed pairs of edge states 
can be explained on the basis of the inability of a state to 
scatter with its time reversed counterpart Q, El • A for- 
mula for the invariant which involves counting the num- 
ber of zeroes of the Pfaffian of the time reversal opearator 
among bands was provided in . This characterization, 
however, does not establish a link between the topology 
in momentum space and the robustness of the edge states 
as observed in numerical studies . In a previous paper 
, we showed that the Zi invariant could be constructed 
as an obstruction in momentum space. The Zi invariant 
was also cast in terms of the more familiar Chern num- 
bers thus leading to a connection between the robust edge 
states and the bulk topology, similar to the one in the in- 
teger quantum Hall effect. Using the additivity of Chern 
numbers, the Zi index for the multi-band case was also 
written down. 

In this paper, we extend the study of Zi invariants to 
three dimensions. We begin with a review of the expla- 
nation of the invariance of the Zi index in two dimen- 
sion elaborating the band touching arguments presented 
in We then extend these methods to study the Zi 



invariants in three dimensions. We find a total of four 
independent Zi invariants, one for each of the faces of 
the Brillouin zone torus and a fourth one which charac- 
terizes the monopole charge in momentum space. A non- 
trivial value of this index implies the existence of a set of 
"trapped monopoles" in momentum space which unlike 
the ordinary case, cannot recombine and are hence sta- 
ble. The monopole charge of an exotic insulator with a 
non trivial fourth Zi index could however change in the 
presence of a magnetic field paving the way for a possi- 
ble experimental observation. We discuss the possibility 
of a three dimensional counterpart of the recently pro- 
posed spin quantum Hall effect We point out that 
while the stability of the edge states in such a model is 
a consequence of the Zi invariants, the stability of the 
quantized response of the spin Hall conductivity is de- 
scribed by a set of three topological invariants in the 
space of twisted boundary conditions, namely, the spin 
Chern numbers [lCJ, which were recently introduced to 
study the stability of the quantized spin Hall response in 
two dimensions. 



TOPOLOGICAL Z2 INVARIANTS 
Z2 invariants in two dimensions 

We begin with a study of the Zi invariant for a ground 
state consisting of time reversed sets of bands in two 
dimensions. First, consider the simplest case when we 
have a single pair of time reversal invariant bands in the 
ground state. Consider a case when the bands do not 
have degeneracies at points other than the ones at which 
such a degeneracy is required by time reversal invariance. 
Note that at degeneracies required by symmetries, in gen- 
eral, the bands can still be identified. The Hamiltonian 
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in this case can be written as 

H = H 1 +H 2 (1) 

where 

H x = [ dkHi{k) 

dk [E 1 (k)\u I {k)){u I {k)\+E 2 (k)\u II (k))(u II (k)\] 
Hi = / dk E n (k)\n(k))(n(k)\ 

n,E„(k)>0 

in which k = (ki,k2),dk = dkidk 2 and E\,E 2 < 0, and 
the Fermi energy, Ef — 0. The time reversal invariance 
of the Hamiltonian implies that 



E 1 (k)=E 2 (-k) 



(2) 



where (— k) = (—ki,—k 2 ) and the bands ui(k),un(k) 
can always be chosen such that 

efluj(fc))) = \u n (-k)), e\u u (k)) = -\ Ul (-k)) (3) 

where O is the time reversal operator. Now imagine that 
we adiabatically change the Hamlitonian H (t) as a func- 
tion of some variable t. As long as the ground state is 
adiabaticaly connected to that of the original Hamilto- 
nian and the bands do not touch, the Hamiltonian can 
always be written as above for some suitably redefined 
bands. The Chern numbers for the bands are then in- 
dividually conserved. The time reversed pair of bands 
have a total Chern number of zero. (This follows from 
homotopy arguments of the form presented in Ij or by 
a direct computation of the Chern number in terms of 
the Berry phase such as by using the TKNN formulae 
[llj . A simple physical argument for this follows from 
the observation that the Hall conductivity of a system 
with time reversal invariance is zero. One can always 
move the bands so that they are the only ones that con- 
tribute to the Hall conductivity. Then using the standard 
result > h follows that the sum of the Chern numbers 
of these bands is zero.) 

Now suppose that a degeneracy occurs at a point 
(to, ko). Due to time reversal invariance, the bands must 
also touch at (to, — fco). It is known that a two level degen- 
eracy is generic while a three or higher level degeneracy 
is not. Further, non-generic three level degeneracies can 
be adiabatically continued to sets of two level degenera- 
cies. Without loss of generality therefore, one can always 
write Hi (t) in the vicinity of the degeneracy points as 



H±(t, k) = m(t,k).a + 



Ei + En 



I 



(4) 



where the two indices of the a matrices correspond to the 
two bands \ui(k)} and \un(k)). Time reversal invariance 
then dictates that 



At the degeneracy point, m(to, k ) — {0, 0, 0}. 

As the parameter t is further changed, the bands may 
split apart again. However, as they separate, the bands 
might now have a new set of Chern numbers whose sum 
has to be the same as before the band collision(zero) [l^ . 
One can visualise the degeneracy points as monopolcs 
which flow in and out of the bands. 

If one sets lo = fn/\m\, then the Chern number ex- 
change between the bands is given by 



n(xo) — ± — / (ujlduj A duj) 
4tt J e . 



(6) 



where Sj is the surface of a small sphere enclosing the 
degeneracy point in the three dimensional space of points 
(t, k x , k y ). From the property of time reversal invariance, 
it follows that the Chern number exchange between a 
band and its time reversed counterpart at the two points 
(to, kg) and (to, — ka) is exactly the same (in magnitude 
and sign). Hence the total Chern number exchanged be- 
tween bands is always an even number. Next consider 
the effects of band collisions for ground states consisting 
of multiple pairs of bands. Two possible kinds of colli- 
sions can occur, the first, a degeneracy between a band 
and its time reversed counterpart, and the second, a de- 
generacy between two bands which are not related by 
time reversal. In the second case, however, the fact that 
the time reversed pairs before and after collision have 
opposite Chern numbers which add up to zero for each 
pair implies that if a band \a) changes its Chern number 
by n a then the Chern number of 0|a) changes by — n n 
These preserve the net Z 2 invariant presented in 



E = | c„|mod2 



(7) 



c„>0 



m(t, —ft) = — fh(t, ft) 



(5) 



The first kind of process also obviously preserves this 
invariant from the previous arguments. Thus we have 
demonstrated that there is an index, E which stays in- 
variant when the Hamiltonian is changed in a continuous 
way such that the Fermi Energy always lies in a gap. 



Zi invariants in 3 dimensions 

Let us parametrize the three torus, T 3 by { — 1 > 
x,y,z < 1}. In three dimensions, the time reversal oper- 
ator, O maps the fiber at the point (x,y,z) to the point 
(—x, —y, —z). Let A be the map on T 3 induced by the 
map O on the fiber bundle of the bands. The planes 
x, y, z = 0, 1 get mapped to themselves and therefore one 
can associate Z 2 invariants with each of them. The planes 
x = y, x — —y, x — z, x — —z, y = z, y = —z are an ad- 
ditional set of planes, which get mapped to themselves. 
We consider separately two different cases : 

1. The bands do not touch at any points in T 3 . 
There are then three independent Chern numbers 
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for each band and the Chern number of the ground 
state is again zero. We can define a set of three 
Z2 invariants E yz , E zx , E xy for each of the planes 
x = 0, y = and z = using the formula for 
the two dimensional case. These are identical to 
the Z2 invariants of the planes x = 1 etc. Thus 
Exy = J2 c my >o c n V \ m od 2 where d^> is the Chern 
number of the nth band in the xy plane and simi- 
larly for E zx and E xy . Now, the Z2 invariants for 
the planes x — y etc. can be obtained as sums of 
the Z2 invariants of the E xy etc, since the Chern 
numbers for these planes can be obtained as sums of 
the Chern numbers of the planes x = 0, y = and 
z = . Let us now perturb H\(x, y, z, t) as before. 
Degeneracies can occur and split into Dirac points. 
These can be thought of as monopoles which can 
recombine after relative displacement through a re- 
ciprocal lattice vector 0] . The Chern numbers for 
the planes are independent and can be changed in- 
dependently. At the same time, because of time 
reversal invariance, the Chern numbers can only 
change in such a way that the Z2 invariants are 
preserved in the process. There are thus only three 
independent Z2 invariants and the Z2 invariants for 
the other planes depend on these. Thus the set of 
all Ground states with TR invariant Hamiltonians 
with no bands crossing the Fermi state and no band 
touching in T 3 can be classified by Z\. 

2. Now let us consider ground states where the bands 
do touch. Let us call the manifolds 

Tl , = [0 > z < 1, -1 > x < 1, -1 > y < 1] (8) 

Tl _ = [-1 > z < 0,-1 > x < 1,-1 > y < 1] (9) 

and similarly for the x and y directions, 3-d half 
tori. Wc further divide these half tori into quar- 
ter tori which are defined by the intersection of 
half tori, such as T 3 = T? f| T 3 • The 

half and quarter tori come in pairs such as {T 3 ^} 
which get mapped onto each other under the pro- 
jection of the time reversal operator. We call these 
"complementary pairs" . For simplicity, we consider 
at first a single pair of time reversed bands and a 
set of generic Dirac like degeneracy points which 
can be thought of as monopoles, which do not lie 
at the surfaces of any of the quarter tori, though the 
arguments can be easily modified to adapt to non- 
generic degeneracies and more complicated cases. 
Time reversal maps these monopoles in quarter tori 
to monopoles in their complementary quarter tori 
such that if a band collides with its time reversed 
counterpart, the charge that flows in at one mono- 
pie to a band is the opposite of the charge that 



flows in at the other monopole. If there is an odd 

number of such monopoles in the half torus T 3 _ , 

2 z 

then the Z2 indices of the plane z = and the 
plane z — 1 arc different, while they are the same 
if there are an even number of monopoles. Since 
these Z2 indices arc homotopically stable, atleast 
one pair of monopoles in the Brillouin zone can- 
not recombine and vanish in the first case. These 
monopoles are thus trapped and hence stable. Fur- 
ther, from the mapping of monopoles to their coun- 
terparts in quarter tori, one can see that if a cer- 
tain half torus has an odd number of monopoles, 
so must all the other half tori. Hence there are 
no additional Z2 invariants. In the case of multi- 
ple bands, one may simply evaluate the monopole 
charge flowing into half the set of filled bands, say 
for instance the set of bands with positive Chern 
numbers, and half the set of those with Chern num- 
ber zero. There are thus four independent Z2 in- 
variants characterizing time reversal independent 
ground states with no bands crossing the fermi sur- 
face. An analysis through K theory seems to sug- 
gest that the fourth Z2 invariant or the monopole 
charge described above, arises from the twisted KR 
theory of S 3 while the remaining three arise from 
the faces of T 3 0. When time reversal symmetry 
is broken, such as in the presence of an external 
magnetic field, the Z2 indices are no longer invari- 
ant, and hence the trapped monopoles may vanish. 
This might lead to an experimental observation of 
such trapped monopoles. 

QUANTUM SPIN HALL EFFECT IN THREE 
DIMENSIONS 

Recently it has been suggested that the spin hall con- 
ductivity can also be quantized in two dimensions and 
various models in which such a quantization occurs have 
since been proposed 0, 0, 0- Here we show that a 
quantized spin Hall response is also possible in three di- 
mensions and discuss the relevance of various topological 
invariants to such an effect. It is well known that in 
two dimensional electron systems with a periodic poten- 
tial, the Hall conductivity in the presence of an external 
magnetic field is quantized to an integer known as the 
Chern number when the Fermi energy lies in a gap of the 
extended states jjjj . The three independent Chern num- 
bers for filled bands in three dimensions play the same 
role in determining the quantized Hall response as the 
single Chern number does in two dimensions 0, 0, 0] . 
The quantized spin Hall current in two dimensions arises 
from equal contributions from the up and down spins, 
J s = 3j(Jt — J j.), each of which has a quantized response 
to the external Electric field, but in opposite directions. 
Thus one expects that in three dimensions, when the 
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Fermi energy lies in a gap, the spin conductivity tensor 
can be written in the form : 



8ir 2 h 



~~ijk 



Gk 



(10) 



where G is a reciprocal lattice vector. When a well de- 
fined Fermi surface exists, the non quantized portion of 
the spin conductivity tensor can be written in terms of 
an integral over the Fermi surface [l4|. A Hamiltonian 
which displays such a quantized response in the k z = 
and k y = planes can be written based on the two dimen- 
sional model presented in || . The stability of the edge 
states is governed by the three dimensional Z 2 invari- 
ants presented earlier, even though the spin Hall conduc- 
tivity itself is not expected to be quantized when terms 
which mix up and down spins are present. An explicit 
formulation of the edge states for the graphene model was 
presented in |20j| and a general connection between bulk 
topology and edge states was discussed in 0. Since a 
non trivial value of the fourth Z2 index is not possible 
without band degeneracies that change the Z2 invariants 
of a pair of bands across a half torus, a quantized spin 
Hall response can only be obtained when this index is 
zero. 

One can instead define three spin Chern numbers fol- 
lowing |ldj for the many body wavefunction of the ground 
state: 



a 



a,0 



4tt 



0/ dip dip 

'08 feet' 



Hi) 



where eijk is the three dimensional antisymmetric tensor 
where a = {1,1} and the variables 6" represent the 
twisted boundary conditions following the notation of 
|lfj|. These are a more appropriate measure of the 
robustness of the spin hall conductivity to the presence 
of spin independent disorder. 

The Z 2 invariant is also discussed in a few other very 
recent preprints. In |2lj . the role of this invariant in a 
one dimensional quantum spin pump is discussed, and 
the role of the index as an obstruction was presented and 
in [22j . the connection between the spin Chern number 
and the Z2 invariant was discussed. 

While this work on 3-d invariants and the 3-d quantum 
spin hall effect was in progress, a preprint [2^] appeared 
which discusses the 2-d and 3-d invariants using results 



from homotopy theory. Our results seem to be in agree- 
ment with the ones obtained in that work using different 
arguments. 
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